data<-read.table("rawvalues_plus2",sep="\t",head=T)
h1<-hist(as.vector(log2(as.matrix(data))),breaks=seq(0,16,0.1))
data<-read.table("rawvalues_plus2_A",sep="\t",head=T)
h2<-hist(as.vector(log2(as.matrix(data))),breaks=seq(0,16,0.1))
data<-read.table("rawvalues_plus2_B",sep="\t",head=T)
h3<-hist(as.vector(log2(as.matrix(data))),breaks=seq(0,16,0.1))
dev.off()

#pdf("histogram_dist_compare_allinone.pdf")
#plot(h1,xlab="Expression after normalization",col="",ylab="Counts",ylim=c(0,10000), xlim =c(0,16),main="EXPR VALUES")
# par(new=TRUE)
#plot(h2,col="",xlab="",ylab="",main="",ylim=c(0,10000))
# par(new=TRUE)
#plot(h3,col="",xlab="",ylab="",main="",ylim=c(0,10000))
#legend("topright", c("PLUS2/PLUS2","HGUA/PLUS2","HGUB/PLUS2","HGUA/HGUA"), title="EXPRESSION VALUES/NORMALIZATION VECTOR",cex=0.8, fill=c("","","",""));
# dev.off()
pdf("histogram_dist_compare_raw.pdf")
plot(h1,xlab="Expression before normalization",ylab="Counts",ylim=c(0,100000), xlim =c(0,16),main="PLUS2/PLUS2")
plot(h2,xlab="Expression before normalization",ylab="Counts",ylim=c(0,100000), xlim =c(0,16),main="HGUA/PLUS2")
plot(h3,xlab="Expression before normalization",ylab="Counts",ylim=c(0,100000), xlim =c(0,16),main="HGUB/PLUS2")
#plot(h4,xlab="Expression after normalization",col="",ylab="Counts",ylim=c(0,10000), xlim =c(0,16),main="HGUA/HGUA") 
dev.off()
